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We report theoretical calculations of high-order harmonic generation (HHG) of Xe with the inclu- 
sion of multi-electron effects and macroscopic propagation of the fundamental and harmonic fields 
in an ionizing medium. By using the time- frequency analysis we show that the reshaping of the 
fundamental laser field is responsible for the continuum structure in the HHG spectra. We further 
suggest a method for obtaining an isolated attosecond pulse (IAP) by using a filter centered on 
axis to select the harmonics in the far field with different divergence. We also discuss the carrier- 
envelope-phase dependence of an IAP and the possibility to optimize the yield of the IAP. With the 
intense few-cycle mid-infrared lasers, this offers a possible method for generating isolated attosecond 
pulses. 
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I. INTRODUCTION 



High-order harmonic generation (HHG), one of the 
most interesting nonlinear phenomena occurring when 
atoms or molecules are exposed to an intense infrared 
laser field, has been widely used for the production of at- 
tosecond pulses in the extreme ultraviolet (XUV) [IJ-Q- 
Due to its great potential for probing ultrafast electronic 
processes, different methods have been used to gener- 
ate isolated attosecond pulses (IAPs). Using carrier- 
envelope-phase (CEP) stabilized few-cycle laser pulses, 
an IAP as short as 80 attoseconds has been generated by 
synthesizing harmonics beyond the cutoff [6J, [7j . Start- 
ing with an elliptically polarized light, the polarization 
gating technique in which HHG emission from the cen- 
tral cycle of the pulse is selected, has also produced an 
isolated 130-as pulse [3Q. In a tight focusing geometry, 
using the so-called spatiotemporal gating, an IAP can be 
generated since different phase-matching conditions can 
be achieved for different ranges of harmonics [IJJ, KLlj . 
Other alternative methods of IAP generation have been 
reported, including confining harmonics generated in a 
narrow temporal window in the leading edge of a laser 
pulse, where a 210-as IAP has been reported [12 |. The 
IAP can also be generate d by op timizing the pressure 
and length of the gas cell [Ia - [l5| or using a spatial fil- 
ter in the far field jlq - KL8l |. Indeed there is a plethora of 
techniques for the production of the IAP, with the idea 
that harmonics be generated from half an optical cycle 
only in a few- or multi-cycle infrared laser pulse. 

Since the harmonic field generated by all atoms or 
molecules within the laser focus co-propagates with the 
fundamental laser field in the medium, as well as possible 
further propagation in the free space depending on the 
experimental setup, the understanding of the observed 
HHG consists of two parts: first, the HHG emission 
from individual atoms (or molecules) through the laser- 
induced dipole; second, the propagation of the fundamen- 



tal and the harmonic fields in the medium and free space. 
There are two fundamental equations to be solved: time- 
dependent Schrodinger equation (TDSE) for the single- 
atom (or single-molecule) response, and Maxwell's wave 
equation for the propagation process. In practice, in the 
first step, the solution of TDSE is normally not employed 
due to computational difficulties. Instead, strong-field 
approximation (SEA) (or Lewenstein model) [191 ] is used. 
However, it is well-known that SEA fails to reproduce 
the observed HHG spectra and thus the reliability of 
the predicted spectra after propagation in the medium 
is questionable. In recent years, a quantitative rescatter- 
ing (QRS) theory has been proposed by our group (20l — 
1221 . which has been corroborated subsequently by others 

23-26]. The single-atom (or single- molecule) response 
obtained from QRS has been shown to be as accurate as 
the one from TDSE, but the calculation is as easy as SFA. 
Taking advantage of this theoretical success, recently we 
have incorporated the QRS theory into the macroscopic 
propagation of harmonic fields in the medium. Until now, 
we have shown that QRS-based macroscopic harmonic 
spectra:_(i) agreed well with TDSE-based HHG spectra 
for Ar [27J; (ii) compared well with experimental HHG 
spectra reported for Ar, N 2 , and CO2 [28H3C| when the 
experimental conditions are well specified. It has been 
demonstrated that this approach can be used to calcu- 
late macroscopic HHG spectra by polyatomic molecules 

31] even though the predictions have not been tested 
against experiments yet. Outside of our group, the QRS 
theory has been applied to study HHG by the two-color 
fields in which the propagation effect is included [32j. 
These applications of the QRS theory are focused on the 
HHG spectra. In this paper, we focus on the analysis 
of attosecond pulse generation which inevitably tests the 
phases of the harmonics obtained in our simulation. 

Recently, Xe has become a favorite candidate for gen- 
erating an intense IAP J33] , studying phase-matching ef- 
fects in the generation of high-energy photons [34|, and 
probing the multi-electron dynamics with high-harmonic 



spectroscopy [35| ■ Ferrari et al. [33| reported the genera- 
tion of a high-energy 160-as IAP using low-order harmon- 
ics of Xe from a CEP-stabilized laser. They used very 
high laser intensity and very dilute gas so that the fun- 
damental field was not severely distorted, but the ground 
state of atom was depleted very quickly in the leading 
edge of the laser pulse. Only low-order harmonics emit- 
ted within one half cycle were used to obtain an IAP. 
Shiner et al. |35| used a 1.8-/im laser with a duration of 
less than two optical cycles to obtain the HHG spectra of 
Xe up to the photon energy of 160 eV. They have shown 
that HHG spectra exhibited strong enhancement above 
about 90 eV. This enhancement is well-known in pho- 
toionization (PI) of Xe due to the presence of a strong 
shape resonance from the 4d shell which, through the 
channel coupling, modifies the partial PI cross section 
of the 5p shell of Xe - a feature attributed to many- 
electron effects. According to QRS, such enhancement is 
anticipated since partial photorecombination (PR) cross 
section (related to photoionization) enters directly in 
the laser-induced dipole. To simulate HHG spectra at 
high-photon energies, multi-electron effects on the laser- 
induced dipoles thus have to be included. Using such 
dipoles in the QRS, we simulate the HHG spectra of Xe 
generated by 1.8-^m lasers by including the macroscopic 
propagation effects. 



In this paper, as a check, we simulate the HHG spectra 
reported in Shiner et al. [351 ] . Mostly we aim at repro- 
ducing the HHG spectra of Xe observed experimentally 
in Ref . [36] , which show nearly continuous photon energy 
distributions (to be called continuum structure) at high 
laser intensities. Such continuum spectra have also been 
observed in molecules, like NO 36]. From our simulation, 
we wish to demonstrate whether IAPs are generated by 
these harmonics. For this, we demonstrate how to se- 
lect different ranges of harmonics to synthesize an IAP 
by using a spatial filter in the far field. This approach is 
different from that in Ferrari et al. [33j, but similar to 
the analysis in Gaarde et al. [13] ■ The rest of this paper 
is arranged as follows. In Sec. II, we briefly summarize 
the propagation equations, wavelet theory for the time- 
frequency analysis, formulas for the attosecond pulse gen- 
eration, and QRS theory including multi-electron effects. 
In Sec. Ill, the theoretical results are presented and an- 
alyzed for different experimental conditions, for IAPs by 
synthesizing harmonic orders from 40 to 80 (H40-H80) 
and H90-H130. We also compare attosecond pulses cal- 
culated using the QRS and the SFA. An analysis of CEP 
dependence of the generated IAP presented at the end of 
this section concludes that it is still possible to obtain an 
IAP even for lasers where the CEP is not stabilized. A 
short summary in Sec. IV concludes this paper. 



II. THEORETICAL METHODS 



Propagation equations of fundamental and 
harmonic fields 



The propagation of the fundamental laser field and 
high harmonics in an ionizing medium has been described 
in detail in Ref. [29(, so we only recall the main equa- 
tions here. The evolution of the fundamental field is 
described by a three-dimensional (3D) Maxwell's wave 
equation |37l.l38j: 



V 2 E x {r,z,t) 



1 d 2 E 1 (r,z,t) dJ ahs {r 7 z 7 t) 

= Mo- 



c 2 dt 2 
-^(l-T^l^r,*,*), 
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(1) 



where E\{r,z,t) is the transverse electric field with cen- 
tral frequency wq. V 2 = V^ + d 2 /dz 2 in cylindrical 
coordinates, where z is the axial propagation direction. 
The effective refractive index is 



r) cS {r, z, t) = ?7o(r, z, t) + r] 2 I(r, z, t) 



Lo 2 (r,z,t) 



2c 2 



(2) 



The linear term rj a = 1 + Si — %f5\ accounts for re- 
fraction (<5i) and absorption (fi\) by the neutral atoms, 
the second term describes the optical Kerr nonlinear- 
ity which depends on the instantaneous laser intensity 
I(t), and the third term contains the plasma frequency 
cj p = [e 2 n c (i)/(e m e )] 1 ' 2 , where m c and e are the mass 
and charge of an electron, respectively, and n e (t) is the 
density of free electrons. The absorption term J^Jt) due 
to the ionization of the medium is given by [J, |39j 



Jabs(t) 



7(t)n c (f)Jpgi(t) 

\Ei(t)\ 2 ' 



(3) 



where j(t) is the ionization rate, and I p is the ionization 
potential. Ionization rates involved in Eq. (j3]) and in free 
electron density n c (t) are calculated using the improved 
Ammosov-Delone-Krainov (ADK) theory [40]. The fun- 
damental laser field is assumed to be Gaussian both in 
space and time at the entrance of the medium, and the 
gas pressure is constant within the medium. 

The 3D propagation equation of the harmonic field is 

'MM 



V 2 E h (r,z,t) 



1 d 2 E h (r,z,t) 
c 2 dt 2 



^o- 



d 2 P{r,z,t) 
~dt 2 : 



(4) 



where P(r,z,t) is the polarization depending on the ap- 
plied fundamental field E\(r, z,t). Here the free-electron 
dispersion is neglected because the frequencies of high 
harmonics are much higher than the plasma frequency. 
In general, the polarization P(r, z, t) is separated into lin- 
ear and nonlinear components, and the linear susceptibil- 
ity x ( w ) includes both linear dispersion and absorption 



effects of the harmonics [43] . The nonlinear polarization 
term P u i(r, z, t) can be expressed as 



P n i(r, z, t) = [n - n c (r, z, t)]D(r, z, t), 



(5) 



where no — n c (r,z,t) gives the density of the remaining 
neutral atoms, and D(r,z,t) is the single-atom induced 
dipole moment. Note that Eqs. (JTJ and Q are solved 
using the Crank-Nicholson routine in the frequency do- 
main. 

Once the harmonic field at the exit face (near field) of 
the medium is computed, the harmonics propagating in 
free space in the far field can be obtained from near-field 
harmonics through a Hankel transformation [44T - I46J . 



In the far field, a spatial filter is used to select the har- 
monics in a prescribed area. In this paper, we assume 
that the filter is circular with a radius r , and is perpen- 
dicular to the propagation direction of harmonics. The 
intensity of an APT or an IAP in the far field is 



ifer(t) 



2irrdr\ / E{(r,u))e iut dcjf 



(10) 



Photorecombination (PR) dipole moment of Xe 
in the QRS theory 



B. Wavelet analysis of attosecond pulses 

A time-frequency representation (TFR) (or spectro- 
gram) of the harmonic field Eh (t) is a simultaneous rep- 
resentation of the temporal and spectral characteristics 
of the harmonics. We perform the time-frequency anal- 
ysis in terms of the wavelet transform of the harmonic 
field StM!: 



A(t,w) 



E h (t')w tyUJ (t')dt' , 



(6) 



with the wavelet kernel w t jM') = y/wW[u)(t' — t)]. We 
choose the Morlet wavelet [471 : 



W{x) = (l/v^e^e 



ix -x jlr 



(7) 



The width of the window function in the wavelet trans- 
form varies as the frequency changes, but the number of 
oscillations (proportional to t) within the window is held 
constant. The dependence of A(t, uS) on the parameter r 
has been tested. The absolute value of A(t,ui) depends 
on r, but the general temporal pattern does not change 
much. In this paper, we choose r = 15 to perform the 
wavelet transform. 

Harmonics emitted at the exit plane (near field) of the 
medium act as a source for the far-field harmonics. In 
order to avoid the complexity of the harmonic spatial 
distribution in the near field (see Fig. 4 in Ref. 29]), 
we calculate A(t,ui) for each radial point in the near field 
and then integrate over the radial coordinate [49| : 

\A aear (t,uj)\ 2 = J 2irrdr\J E h (r,t')w t ^(t')dt'\ 2 . (8) 

To demonstrate the divergence of harmonics, we preform 
the TFR for each radial point in the far field. 

The spectral filter used to select a range of harmonics 
(wi - U2) could affect the generation of attosecond pulse 
trains (APTs) or IAPs. Theoretically we can obtain the 
total intensity of an APT or an IAP in the near field as 
following [51j : 



-Micar \t) 



The single-atom induced dipole moment D(t) in Eq. 
([5]) is obtained by the QRS theory. It can be expressed 
in energy (or frequency) domain as following [22], |52j : 



D(u) = W(u)d(u) 



(11) 



2irrdr\ / E h (r,uj)e lu ' t duj\' 2 



(9) 



where d(u) is the PR transition dipole moment and W{uj) 
is the microscopic wave packet. In QRS theory, W{u) is 
determined by the laser field and can be accurately calcu- 
lated from SFA, and d(ui) is the transition dipole between 
the initial and final states of PR (or PI) . When the multi- 
electron effect is not important, the transition dipole can 
be calculated using the single-active electron (SAE) ap- 
proximation. However, the transition dipole is easily 
generalized to include many-electron effects, as routinely 
done in PI theory of atoms and molecules. Thus to in- 
clude many-electron effects in d(u)), multi-channel cal- 
culations such as many-body perturbation theory, close- 
coupling method, R-matrix method, random-phase ap- 
proximation, and many others can all be employed for 
such purpose. Since PI of Xe has been well studied, we 
obtain d(u>) semi-empirically. The major many-body ef- 
fect for PI of Xe from 5p shell occurs at photon energy 
where 4d shell is open. Thus below about 60 eV, the 
transition dipole from 5p can be obtained from a single- 
electron model. This gives the magnitude and phase of 
the transition dipole. At higher energies, effects from 
the 4d shell on the transition dipole of 5p becomes im- 
portant since PI cross section of Xe from 4d has a large 
and broad shape resonance around 100 eV. The intershell 
coupling will enhance the d(u)) for 5p near and above 90 
eV. Such enhancement has been calculated in Kutzner 
et al. [531 ] using the relativistic random-phase approx- 
imation (RRPA). In our calculation, the phase of d(ui) 
is taken from the 5p shell under the SAE approxima- 
tion [29| while the magnitude is taken from Ref. [53j . 
This approximation does not change the temporal struc- 
ture of attosecond pulses (shown later) since the phase 
of D(u) is dominated by the phase of the wave packet 
W{ui), We comment that in QRS the induced dipole is 
given in the energy domain, thus the calculation is simi- 
lar to the time independent theory used in PI which has 
been well-established in the last 30 years. 
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FIG. 1: (Color online) Macroscopic HHG spectra of Xe in an 
1825-nm laser, for (a) CEP=0 and (b) CEP averaged. Laser 
intensities are indicated in units of Io=10 14 W/cm 2 . See text 
for additional laser parameters and the experimental arrange- 
ment. 



III. RESULTS AND DISCUSSION 



A. Macroscopic HHG spectra of Xe at different 
laser intensities 



HHG spectra of Xe extended to the photon energy of 
over one hundred electron volts using 1.8-/xm lasers with 
the pulse duration of few optical cycles have been re- 
ported recently |3a, |36| . 

In Fig. [TJ we show the calculated HHG spectra of Xe 
exposed to a 14-fs (FWHM), 1825-nm laser. The laser 
beam waist is 100 fim. A 1-mm-long gas jet with the 
pressure of 30 Torr is placed at the laser focus. The har- 
monics are detected after a slit with a width of 190 /im 
and placed 455 mm behind the focus. These parame- 
ters are chosen to be close to those in the experiment of 
Trallero-Herrero et al. 36]. For the present purpose we 
analyze HHG spectra obtained from our theoretical simu- 
lations at two laser peak intensities 0.5 x 10 14 W/cm 2 and 
l.OxlO 14 W/cm 2 , which are below and above the critical 
intensity for Xe at ~ 0.87x 10 14 W/cm 2 [|oj], respectively. 
Here the critical intensity is defined with respect to the 
static electric field where an electron can escape over the 
top of the field-induced potential barrier classically. 

We show the macroscopic HHG spectra for CEP=0 
in Fig. [Ha). The two laser intensities present different 
characteristics of harmonics. For the low intensity, the 



harmonics are very sharp, i.e., the valley between the 
neighboring odd harmonics is very deep. At high inten- 
sity, the valley is very shallow, i.e., the spectrum shows 
a continuum structure. Furthermore, the harmonics are 
not exactly at odd orders due to the blue shift of the fun- 
damental field. Note that the spectrum rises above about 
H90 is due to the intershell or many-electron effects dis- 
cussed in Sec. II C. Since a few-cycle laser pulse is ap- 
plied, the HHG spectra have a strong CEP dependence. 
In Fig. QTb) , we show the CEP averaged HHG spectra. 
The main characteristics of harmonics remain the same 
except that the harmonic spectra are much smoother. 
The CEP is fixed at zero in the following sections unless 
otherwise stated. 



B. Spatiotemporal evolution of the fundamental 
laser field 

To understand the different spectral features in Fig. [TJ 
we inspect the fundamental field in the ionizing medium. 
The spatiotemporal intensity profile and on-axis electric 
fields of the laser pulse at the entrance and the exit of the 
gas jet are shown in Fig. [3] The laser peak intensity is 
l.OxlO 14 W/cm 2 , which would give an ionization prob- 
ability of ^35% at the end of laser pulse for Xe accord- 
ing to an empirical ADK formula in barrier-suppression 
regime [40]. While the electric field at the entrance has 
a good Gaussian shape both in time and space, it is 
strongly reshaped during the propagation in the ionizing 
medium. At the exit it shows positive chirp in time (blue 
shift in frequency) [see Fig. Hfc)] and defocusing in space 
[see Fig. [5[b)]. We have also checked the fundamental 
field with laser peak intensity of 0.5x 10 14 W/cm 2 . It al- 
ways maintains Gaussian spatial distribution and there 
is no blue shift because the ionization probability is very 
low. The reshaping of the fundamental field at high in- 
tensity is responsible for the continuum structure in the 
HHG spectra in Fig. [TJ Note that similar results have 
been obtained by Gaarde et al. [18| using a 750-nm laser 
interacting with Ne gas. 



Time-frequency analysis of harmonics in the 
near and far fields 



Harmonic generation is a temporal coherent process 
and can be better understood if we study it in terms 
of its emission time. In this subsection we examine the 
time- frequency representation (TFR) of harmonics in the 
near and far fields for the low and high laser intensities. 

For each harmonic order q, it is known that the phase 
can be expressed as [54 ]: 



$l(r,z,t) = -a$I(r,z,t), 



(12) 



where I(r,z, t) is the spatiotemporal intensity of the fun- 
damental laser field. The proportional constant on—s, l 
depends on "short" (S) or "long" (L) trajectories. The 
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FIG. 2: (Color online) Spatiotemporal intensity profile of the laser pulse at (a) the entrance and (b) the exit of Xe gas jet. Laser 
intensity at the focus is 1.0 xlO 14 W/cm 2 (assumed in the vacuum) and CEP=0. (c) Evolution of the on-axis electric field at 
the entrance (solid line) and the exit (dot-dashed line). The laser field becomes chirped during the propagation. For sub-cycle 
dynamics analysis, we use the label B t , with t=-l, -0.5, 0, and 0.5 (in units of optical cycles) to indicate the approximate 
half-cycle where the electron is born. Note that t is defined within the half cycle only. 



phase can also be expressed in terms of the pondero- 
motive energy U„ an d the electron excursion time r/ : 
<J>* cs —fiiUpT? [53| , where the coefficient /% for the 
"short" trajectory is much smaller than for the "long" 
trajectory. The electron excursion times for the two tra- 
jectories are rj ~ T/2 and t\ w T (T is the laser pe- 
riod) [56]. It shows that the phase grows with the cu- 
bic power of the wavelength. The curvature of the phase 
front caused by the radial variation d^f{r)/dr makes the 
harmonic beam divergent. The divergence of "short"- or 
"long" -trajectory harmonic is determined by either Aa- 
or AI(r). 



1. Harmonics in the near field 

The TFR, \A ncal (t, w)| 2 , calculated from Eq. flS), are 
shown in Figs. |3t[a) and (d) for harmonics above H40 
at two laser intensities, collected at the exit face of the 



gas jet (near field). In Fig. [3Ja), the symbols S and L 
are used to indicate the first (earliest) group of harmon- 
ics generated. Here S (L) stands for "short" ("long")- 
trajectory harmonics that have positive (negative) chirp. 
These harmonics are from electrons born at t=-l (in units 
of optical cycles), i.e., B_i to indicate born time at t=- 
1, in the leading edge of the pulse [see Fig. HJc)]. In 
the following, the electron born time t (in units of op- 
tical cycles) is indicated by B t in the figure, while the 
harmonic emission time is read off from the horizontal 
axis of the figure, one for the "short", and the other 
for the "long" trajectory. In this paper the time is al- 
ways defined in moving coordinate frame [29[. At the 
low intensity in Fig. G^a), we can see that both S and L 
contribute to harmonics generated from electrons born at 
t=-l, -0.5, 0, and 0.5. In other words, harmonics are gen- 
erated by electrons born over four half cycles. Note that 
Tate et al. [571 ] have shown that harmonics generated 
by mid-infrared lasers had large contributions from elec- 
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FIG. 3: (Color online) Top row: Time-frequency representation (TFR) of harmonics in the near field. Middle row: TFR for 
on-axis (r=0 mm, divergence: mrad) harmonics in the far field. Bottom row: TFR for off-axis (r=l mm, divergence: 2.2 
mrad) harmonics in the far field. Far-field position is at z=455 mm, and laser intensity (CEP=0) along each column is indicated. 
Electrons are released at each half cycle, labeled by Bt, with t=-f, -0.5, 0, and 0.5 as in Fig. [2] For each B t , electrons can 
follow a "short" (S) or "long" (L) trajectory to recombine with the ion to emit harmonics. For each harmonic, the emission 
time can be read from the time axis. For each Bt, the emission time for each off-axis harmonic is delayed with respect to the 
corresponding on-axis harmonic, e.g., compare (b) vs (c), and (e) vs (f). All the TFRs have been normalized. 



tron trajectories even longer than the "long" trajectories 
in single-atom response, which has also been confirmed 
in our calculation (not shown). But these trajectories 
are all eliminated during the propagation in the medium 
since their phases are very large. For low intensity, the 
propagation in the medium cannot eliminate contribu- 
tions from "long" trajectories. 

The same TFR analysis for the high intensity is shown 
in Fig. [3jd). Higher harmonic cutoff from each burst is 
easily seen since the intensity is twice higher. Comparing 
to Fig. |3{a) , there are no contributions to the harmonics 
from the "long" trajectories for electrons born at t=-l 



and -0.5, i.e., from the leading edge of the pulse. Since 
the laser intensity is twice higher, the phase of each har- 
monic is also twice higher (also see Fig. 17 in Ref. |58j 
and Fig. 1(A) in Ref. j59|]), thus resulting in cancela- 
tion of contributions from the "long" trajectories. For 
electrons born at the falling edge of the pulse, due to 
the blue shift (thus shorter wavelength) and reshaping 
(thus lower intensity) the phases of harmonics due to the 
"long" trajectories are smaller and they can survive after 
propagation in the medium, for example, for electrons 
born at t=0, and 0.5, see Fig. [3Jd). 



2. On-axis harmonics in the far field 



D. Spectral and spatial filtering in the generation 
of attosecond pulses 



In Fig. [3(b), the TFR is shown for r=0 mm in the far 
field (455 mm after the laser focus) . At low intensity, the 
emission from "short" trajectories born at different times 
have the similar small divergence, and after propagation 
in free space they all survive along the axis in the far field. 
Interference between "short" -trajectory harmonics from 
each half cycle leads to enhancement in odd harmonics 
and suppression in even harmonics (see Fig. 18 in Ref. 
|60| ) , and resulting in a big contrast between an odd har- 
monic and neighboring harmonics shown in Fig. QJa) f° r 
the spectra obtained with a slit. At high intensity, only 
harmonics from the "short" -trajectory electrons born at 
t=-l survive (the next one at t=-0.5 is much weaker), see 
Fig-EIe). This would result in a nearly continuum spec- 
tra and a potential for generating an isolated attosecond 
pulse. 



-axis harmonics in the far field 



In Figs. [3(c) and (f), the TFR is shown for r=l mm (di- 
vergence: 2.2 mrad) in the far field. Each off-axis burst 
has an obvious time delay with respect to the on-axis 
burst because it travels a longer distance in free space. 
At low intensity, harmonics from "long" trajectories ap- 
pear on each burst since they have large divergence [see 
Fig. [3(c)]. At high intensity, the "short" trajectories 
contribute to bursts B_o.s and Bo [see Fig. [3(f)]. They 
appear to come from the pulse reshaping, see Fig. [2(b) 
showing laser peak intensity shifting to region away from 
the propagation axis. They experience larger A/(r) with 
respect to "short" -trajectory electrons born at B_i at the 
leading edge. Fig. [3(f) shows that a continuum spectra 
from a "short" trajectory is generated for electrons born 
at t=-0.5. 

Note that attochirp (emission time varying with har- 
monic order) [59|,[6l| of "short"- or "long" -trajectory har- 
monics exists even after propagation. They may be com- 
pensated using a "plasma compressor" [59j because free 
electrons induce a negative group velocity dispersion, or 
by thin filters with linear negative group velocity dis- 
persion [62j. But attochirp is inversely proportional to 
laser wavelength [631 ]. This implies that one can select a 
broad range of harmonics to synthesize a short attosec- 
ond pulse using an 1825-nm laser (will be shown next). 
The harmonic emission of "short" trajectory in the far 
field in Figs. [3(e) and (f) varies with time or radial dis- 
tance. This provides the possibilities to generate IAPs 
using the different ranges of harmonics on or off axis. 
We will only show the spectral and spatial filters applied 
on axis in the far field in the following. 



A spectral filter is usually used to synthesize attosec- 
ond pulses. In this section we also study how the at- 
tosecond pulses are manipulated through spatial filtering. 
Fig. @}a) displays the intensity profile of an XUV light 
by synthesizing H40-H80 at the near field generated by 
laser intensity of 0.5xl0 14 W/cm 2 . The intensity of the 
attosecond pulses I nea r(t) is calculated by using Eq. ©. 
The time-frequency analysis of these harmonics has been 
given in Fig. [3(a). Besides attosecond bursts occurring at 
each half optical cycles, which can be attributed to har- 
monics resulting from "short" trajectories, we observe 
other pulses in between which are attributed to contri- 
butions from "long" trajectories. The main peaks from 
the "short" trajectories are labeled by 1, 3, 5, and 7 in 
the figure, while those in between (2, 4, 6, and 8 are not 
labeled) arc from "long" trajectories. The attosecond 
pulses thus generated show a poor periodicity in time, 
see Fig. E(a). 

If the XUV light is synthesized at the far field, in par- 
ticular, by introducing a spatial filter, then it may be 
possible to remove harmonics resulting from the "long" 
trajectories. In Fig. E(b), the intensity distributions of 
the synthesized light in space in the far field are shown. 
They are obtained from the near-field harmonics by fur- 
ther propagation in free space. The peaks 2, 4, 6, and 8 
[not shown in Fig. 0(a) explicitly] are attributed to "long" 
trajectories. They are indicated in Fig. [4(b) showing that 
they are distributed far from the propagation axis. By 
using a spatial filter (indicated by a solid line in red, with 
a radius ro=100 /im) to select harmonics generated near 
the axis only, as shown in Fig. H(c) by using Eq. (fTU|) to 
calculate I/ ar (t), well-behaved APTs are then obtained. 
We comment that the time delay between off-axis and 
on-axis harmonics leads to the curved spatial distribu- 
tion in Fig. [4(b), and it can be understood mathemati- 
cally since each harmonic behaves like a Gaussian beam, 
and the geometric phase of each harmonic is proportional 
to r 2 along the transverse direction (see Fig. 4 in Ref. 
[23]). The traveling distance of off-axis harmonics can be 
compensated using a reflecting mirror to refocus the har- 
monic beam or by a detector with a curved surface. In 
principle, this compensation becomes important to re- 
duce the duration of attosecond pulses when a spatial 
filter with a large radius is applied. In this paper, the 
radius of the spatial filter is chosen to be small enough 
to avoid this curvature effect. 

Next we use the same range of harmonics (H40-H80) 
generated by the laser intensity of l.OxlO 14 W/cm 2 to 
synthesize attosecond pulses in the near field. Referring 
to Fig. [3(b), the "short" trajectories dominate the har- 
monic generation in the leading edge of the laser, while 
the "long" trajectories dominate the harmonic generation 
in the falling edge. The synthesized XUV light, shown 
in Fig. 2(d) indeed reflects this point where the first two 
peaks occur at multiples of half optical cycles, while the 
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FIG. 4: (Color online) First column: Intensity (or envelope) of attosecond pulses in the near field, synthesized from the 
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trajectories) have large divergence, or at large r. The odd bursts (not labeled) have smaller divergence. There is a time delay 
between off-axis attosecond pulses compared to on-axis ones. Last column: Intensity of attosecond pulses in the far field using 
a spatial filter with a radius ro=100 fjxa (shown by the solid line in red in each middle-column frame). 



last four peaks are not. In Fig. life), the spatial distribu- 
tion of the synthesized XUV light in the far field indeed 
supports this description. By using a spatial filter (indi- 
cated by a solid line in red, with a radius ro=100 /an) to 
select only "short" trajectories, as shown in Fig- Sff), a 
nice IAP with a duration of 270 as is obtained, accompa- 
nied by a weak sub-pulse with a much weaker intensity. 
This demonstrates the generation of IAPs using spatial 
filtering. A similar mechanism of IAP generation has 
been proposed by Strelkov et al. [6J, |65[ using the har- 
monics in the plateau region generated by the Ar gas 
with very high pressure. 

The TFR in Fig. [3(e) shows considerable on-axis emis- 
sion above H80 at burst B_o.5- We use H90-H130 to 
generate attosecond pulses in the near field in Fig. 3(g) . 
Both bursts have considerable contributions from "short" 
trajectories. In the far field [see Fig. 3(h)], they show dif- 
ferent divergences as discussed before. Finally, we obtain 
an IAP with a duration of about 170 as in Fig. 3(i) with 
a spatial filter. The intensity of the IAP is about 1/8 as 



that in Fig. 3(f) due, not only to the larger divergence of 
"short" -trajectory harmonics born at B_o.s than at B_i, 
but also the lower harmonic intensity of H90-H130 than 
that of H40-H80. On the other hand, the duration of the 
IAP is decreased. Similar mechanism of IAP generation 
has been proposed by Gaarde et al. [13, LL8| using har- 
monics in the cutoff region by a 750-nm laser exposed on 
Ne gas. 

We next check how the IAP generation works at higher 
intensities, say at 2.0xl0 14 W/cm 2 . Fig. [5] should be 
compared to Figs. 3(d)-(i) directly. For the synthesized 
H40-H80 pulse, Fig. [5(a) shows that the harmonics are 
emitted about half an optical cycle earlier than the one 
at half the intensity (l.OxlO 14 W/cm 2 ). Fig. 13(b) shows 
that the pulses generated at the falling edge of the laser 
pulse have large divergence and thus they tend to come 
from "long" trajectories. In fact, this portion of the pulse 
does not have good periodic time dependence. Fig. [5(b) 
also shows that only the pulse emitted at t=-l (in units 
of optical cycles) is near the axis, thus a filter selecting 
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FIG. 5: (Color online) Top row: Attosecond pulses synthesized from H40-H80, (a) at the near field, (b) its spatial distribution 
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IAP for H40-80 and for H90-130 are indicated by B t for t=-1.5 and -1, respectively. Laser intensity is 2.0xl0 14 W/cm 2 and 
CEP=0. 



near axis harmonics results in an IAP, as illustrated in 
Fig. EJc). The IAP has a duration of 260 as. For pulses 
synthesized from H90-H130, Fig. EJd) shows that there 
are two bursts emitted at t=-0.5 and 0, and their lateral 
profiles in the far field are shown in Fig. Ete) . By using a 
filter, an attosecond pulse of 150 as can be obtained. In 
this case, the IAP intensity does not increase since the 
fundamental laser field is much reshaped and the har- 
monic has a much bigger divergence in the far field in 
comparison with 1.0 xlO 14 W/cm 2 . It is concluded that 
ionization gating still works at a higher laser intensity, 
and it is more efficient to select bursts that are born in 
the leading edge before the laser field starts to be de- 
pleted and blue shifted. 



Another question arises is whether the strength of at- 
tosecond pulses can be improved by increasing gas pres- 
sure. For "weak" fields and low pressure, Shiner et al. 
[66j have shown experimentally that the harmonic yield 
increases quadratically with the pressure. For "high" 
fields addressed here, laser pulse reshaping is important, 
we have confirmed theoretically that increasing the gas 
pressure while maintaining the same intensity would not 
always increase the harmonic yields of Ar [67|. On the 
other hand, the pressure effect on the high harmonics of 
Xe at "high" field has been studied in Ref. [36|]. But its 
effect on the attosecond pulse wasn't examined in this 
paper. 



E. Far-field position dependence of isolated 
attosecond pulses 

The position of the spatial filter in the far field can 
be easily adjusted in an experiment. Here we show the 
change of attosecond pulses with the far-field position. In 
Fig. [51 attosecond pulses synthesized (H40-H80) at two 
other positions z=100 mm and 900 mm are given. This is 
to be compared with the ones at z=455 mm in Figs.rjje) 
and (f). Using the same filter (indicated by the solid 
line in red for fixed ro=100 yum), the attosecond pulses 
generated are shown in Figs. HJb) and (d). The width 
of the main burst does not change much with z, but the 
strength of the satellite peak is reduced. Of course this is 
achieved at the expense of decreasing the strength of the 
attosecond pulse. This has also been shown by Gaarde 
and Schafer [17J where an IAP was selected by moving 
the reflecting mirror further from the laser focus. 



F. Comparison between QRS and SFA in modeling 
propagation effect 

In the last two decades, the strong-field approximation 
(SFA), which is in the frame of the SAE approximation, 
has been widely used to predict the temporal structure 
of attosecond pulses even though SFA is unable to ex- 
plain the observed harmonic spectra precisely in general. 
In the present calculation, we use QRS in the propa- 
gation calculation. For single-atom response, QRS has 
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been tested against TDSE, both for the magnitude and 
phase, as documented in Le et al. [52i |. for example. In 
QRS, the wave packet is obtained from SFA, including 
the phase. The transition dipole d(ui) in Eq. ([TTt . also 
introduces a phase. In SFA, this phase is a constant, ei- 
ther real or pure imaginary (depending on the symmetry 
of the ground state) and independent of the harmonic 
order. In QRS, the transition dipole moment is a com- 
plex number in general. From PI theory, however, it is 
known that the phase of the transition dipole does not 
change much with the photon energy. Thus the phases 
of the harmonics calculated from QRS and SFA do not 
differ significantly. Since the phases of the harmonics are 
much more important in synthesizing attosecond pulses 
[68| , this explains why propagation theory based on SFA 
has been so successful in explaining the generation of at- 
tosecond pulses, in spite of its failure in predicting or 
explaining the observed harmonic spectra. In this sub- 
section, we support this analysis with actual results from 
simulations. 

In Fig. Eta) the HHG spectra of Xe obtained from SFA 
(within the SAE approximation) and QRS (including 
multi-electron effects) using the laser parameters given in 
the captions are shown. Clearly the spectra differ greatly. 
In Figs. Etb) and (c) the synthesized (H40-H80) attosec- 
ond pulses at the near field and the far field are shown. 
Clearly the results from the two calculations are essen- 
tially identical (after an overall normalization), in spite 
of the large differences in the harmonic spectra. We have 
checked some other cases and found that the temporal 
structures of the attosecond pulses from the two theories 
are always very similar. Larger differences than those 
shown in Figs. EJb) and (c) are expected if a wider range 
of harmonics are used or if the spectra from the two the- 
ories differ much more, but the general conclusion is cor- 
rect. 



G. CEP dependence of isolated attosecond pulses 

The selection of an IAP by a spatial filter in the far 
field discussed above is only for a single CEP, and thus 
only useful if the laser is CEP-stabilized (has not been 
achieved for 1.8- /mi lasers yet). To check if the method 
can be used for lasers that are not CEP-stabilized, we 
investigate the CEP dependence of the IAP generation. 

In Figs. EJa) and (b), we show the contrast ratio be- 
tween the intensities of the strongest satellite and the 
strongest attosecond burst, and in Figs. [SJc) and (d), 
we show the peak intensity of the strongest attosecond 
burst, as the CEP is varied, for the two laser intensities 
indicated. A good IAP is to have high peak intensity for 
the main peak and weak satellites. From Figs. H£c) and 
(d), we note that at the CEP's where the strongest at- 
tosecond bursts have high peak values, the contrast ratios 
shown in Figs. Ufa) and (b) at these CEP's are always 
small. In the meanwhile, when the contrast ratio is large, 
the strongest attosecond burst is always weak. Thus it is 



possible to generate single attosecond pulses even when 
the CEP of the driving laser is not stabilized. This ex- 
plains the success why the first single attosecond pulses 
were generated using few-cycle laser pulses that were not 
phase-stabilized (16j |. 



IV. CONCLUSIONS 



In this paper, we have studied the generation of 
isolated attosecond pulses (IAPs) using few-cycle mid- 
infrared lasers at large intensities near and above the 
critical intensity of Xe. The calculations are based on the 
QRS theory where many-electron effects are included in 
the single-atom induced dipole moment; specifically, by 
including the coupling of the inner 4d shell of Xe in the 
partial 5p photorecombination transition dipole matrix 
element. The effect of the medium on the fundamental 
and harmonic fields is obtained by solving the Maxwell's 
wave equations. The modification (or reshaping in space 
and time) of the fundamental field is due to its nonlinear 
interaction with the medium includes dispersion, plasma 
effect, and Kerr nonlinearity. We have investigated the 
spatiotemporal evolution of the fundamental laser field in 
detail, and found that its reshaping is responsible for the 
continuum structure in the HHG spectra. This conclu- 
sion is carried out in terms of the time-frequency analysis 
of harmonics in the near and far fields. 

Since the divergence of harmonic emission from differ- 
ent half cycles is varied due to the blue shift and defo- 
cusing of the fundamental laser pulse (or complicated re- 
shaping), we have shown that isolated attosecond pulses 
can be generated by synthesizing H40-H80 or H90-H130, 
selected by a spatial filter centered on the propagation 
axis in the far field. The mechanism of IAP generation in 
this paper could be called as "ionization gating" . It works 
for a loosely focused laser at high laser intensity (above 
the critical intensity) , which is reshaped as it propagates 
through the medium with a moderate gas pressure. A 
similar approach has been discussed by Gaarde et al. 
p7l [l8[ using a 750-nm laser interacting with 135-Torr 
Ne gas. We have found that it is easier to reshape the 
fundamental field using a long-wavelength laser with a 
moderate gas pressure (~ 30 Torr). The extended har- 
monic cutoff of Xe leads to a broad range of harmonics 
available for IAP generation. This approach is also dif- 
ferent from Ferrari et al. [33j where low harmonics (~ 
30 eV, which is equivalent to H40 in this paper) are used 
to generate the IAP. In addition, we have discussed the 
possibilities of improving the intensity of the IAP, such as 
changing far-field position, increasing laser intensity, and 
increasing gas pressure. We have shown that the method 
is very robust and an IAP can be generated even if the 
laser CEP is not stabilized. 
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